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Abstract 

We present a search for the rare decay B° -> K*° vu, using a data sample of 492 fb 1 collected 
with the Belle detector at the KEKB e + e~ collider. Signal candidates are required to have an 
accompanying B meson fully reconstructed in one of the hadronic modes and signal-side particles 
consistent with a single K*° meson. No significant signal is observed in the data sample and we 
set a 90% confidence level upper limit of 3.4 x 1CT 4 on the branching fraction. 

PACS numbers: 13.25.Hw, 14.40.Nd 
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The decays B — > K^vv and B — > K^l + l~ proceed through the flavor-changing neutral- 
current processes, which are sensitive to physics beyond the Standard Model (SM) associated 
penguin loops. In the SM, the dominant diagram of these decays is the penguin process 
shown in Figure ^ a). The SM branching fractions are estimated to be around the 10 _5 /10 -6 
level for the K* / K modes based on next-to-leading-order calculations pj. Calculation of the 
decay amplitudes for B — > K^uV is particularly clean theoretically, owing to the absence 
of long-distance interactions that affect the charged-lepton channels B — > K^l + l~ . New 
physics such as SUSY particles or the effect of a possible fourth generation could potentially 
contribute to the penguin loop or box diagram (Figure H b)) to enhance the branching 
fractions Q. Reference also discusses the possibility of discovering light dark matter in 
b — > s transitions with large missing momentum. 



b 
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FIG. 1: The quark-level diagrams for B — ► K*uv decays. 

Experimental measurements of the b — > s transitions with two charged leptons are 
in good agreement with SM calculations []]]. Further investigation of the forward-backward 
asymmetry in B — > K*l + l~ 0| prefers the SM prediction although the statistics are still 
limited. Due to the challenge of cleanly detecting rare modes with two final-state neutrinos, 
only a few studies of K^vV have been carried out to date |5[. In this paper, we report our 
first search for the decay B° — > K*°uV using a data sample of 492 fb" 1 integrated luminosity 
recorded at the T(4S f ) resonance, corresponding to 535 x 10 6 B- meson pairs. 

The Belle detector is a large-solid-angle magnetic spectrometer located at the KEKB 
collider and consists of a silicon vertex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cherenkov counters (ACC), a barrel-like arrangement 
of time-of-flight scintillation counters (TOF), and an electromagnetic calorimeter comprised 
of CsI(Tl) crystals (ECL) located inside a superconducting solenoid that provides a 1.5 T 
magnetic field. An iron flux-return located outside of the coil is instrumented to detect K\ 
mesons and to identify muons (KLM). The detector is described in detail elsewhere Q. 

One of the B mesons in the event is fully reconstructed as the tag-side B candidate (-B tag ). 
The rest of the particles are assumed to be products of the signal-side B meson (B sig ). The 
B Ug candidates are reconstructed in one of the following decays: B° — > D^~n + , D^*'p + , 

D^-af, and D^-D { * )+ . The D~ mesons are reconstructed as D — > K^tt , K^n n°, 
Kg7T~7T + 7T~ , K + 7t~tt~, and K + 7t~7t~7t° . The D*~ mesons are reconstructed as D°ir~, and 
the following decay channels are included for the D° mesons: D° —> K + n~, K + n~n°, 
K+n-n+n-, K° s ir°, K° s n-7i + , K s 7r-n + n and R- K + . Furthermore, the D+ -> K° S K + and 
K + K~tt + decays are used for mesons and D* + — > decay are also selected. A 

candidate B tag meson is selected with the beam constrained mass M^ c = ■sj E^ eam — p 2 B and 
the energy difference AE = E B — -Ebcam, where E B and p B are the reconstructed energy 
and momentum of the _B tag candidate in the T(45 r ) center-of-mass (CM) frame, and E heam 
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is the beam energy in the same system. The requirements on candidate B ta , g mesons are 
M bc > 5.27 GeV/c 2 and -80 MeV < AE < 60 MeV. If there are multiple B tag candidates 
in one event, the candidate with the smallest \ 2 based on the deviations from the nominal 
values of AE, the D meson mass, and the mass difference between the D* and the D (for 
the candidate with a D* in the reconstruction) is chosen. 

Once a B tag candidate is reconstructed, the remaining particles are used to reconstruct 
a B° — > K*°vV event. The K*° candidate is reconstructed by two charged tracks with 
opposite charge, and one of the tracks should have a kaon likelihood greater than 0.6. The 
kaon likelihood is defined by Tlx = £k/(£k + Ar)> where Ck denotes a combined 
likelihood measurement from the ACC, the TOF, and a dE/dx from the CDC for the K ± 
(tt^) tracks. The daughter K ± and tt t are required to have a maximum distance to the 
interaction point (IP) of 2 cm in the beam direction (z) and of 5 cm in the transverse plane 
(r-(j)). The invariant mass of the K*° candidate should be within a ±75 MeV/c 2 window of 
the nominal K*° mass. No other charged tracks or tt° candidates are allowed in the event, 
while a pair of photons with an invariant mass within ±18.5 MeV/c 2 of the nominal it 
mass where each photon has an energy greater than 50 MeV is considered a tt° candidate. 
A B S i g candidate is selected according to the variable -Eecl = E to t — -^rec, where E tot and 
E rec are the total visible energy measured by the ECL detector and the measured energy 
of reconstructed objects including the B tSbg and the signal side K*° candidate, respectively. 
A minimum threshold of 50 (100,150) MeV on the cluster energy is applied for the barrel 
(forward endcap, backward endcap) region of the ECL detector. The signal region is defined 
by E E cl < 0.3 GeV and the sideband region is given by 0.45 GeV < E EC l < 1-2 GeV. 

The dominant background source is generic B°B° decays. As shown in Figure El a 
Fisher discriminant containing three input variables (-P£., -Pmi SS , and M^ iss ) is introduced 
to suppress the background, where P£«, -P^iss^ and M^ iss are the momentum of the K* 
candidate in the CM frame, the missing momentum in the CM frame, and the squared 
missing mass. The missing momentum and missing mass are calculated using the momenta of 
the reconstructed B tag and K*° candidate. Based on a figure of merit study, the events with 
a Fisher discriminant value greater than —4.9 are rejected from the analysis. Furthermore, 
the cosine of the angle between the missing momentum in the laboratory frame and the beam 
pipe direction is required to be within the range —0.86 and 0.95. These criteria suppress 
events with particles produced along the beam pipe. The contributions from continuum 
background e + e~ qq (q = u, d, c, s) and other rare B decays such as B° — > K* Qr y are 
expected to be small. Based on Monte Carlo (MC) simulations, the selection efficiency 
on the signal side is estimated to be 11%, while the B tagging efficiency is 0.087%. The 
reconstruction efficiency in this analysis is equal to the product of the selection efficiency 
and the B tagging efficiency 

The contribution from each background source is examined with large MC samples. The 
expected yields of background events in the signal region and sideband region are 4.8±1.5 and 
19 ± 3 events, respectively. Further details are listed in Table HI From the full experimental 
sample 13 events fall in the signal region and 19 events fall in the sideband region. 

The signal yield is extracted by a fit to the -Eecl distribution. The signal probability 
density function (PDF) is modeled by a smooth histogram obtained from a signal MC 
sample, and a second order polynomial is used to describe the background distribution. An 
extended likelihood function is introduced 

-(N S +N B ) N 

£ = m II [NsPs(E* ECh ) + N N P B (E ECL )} , (1) 

i=i 
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where Ps and Pb denote the signal and background PDF. The signal yield, background yield, 
and total number of events in the fit are given by Ns, Nb, and N, respectively. We obtain 
a signal yield of 4. 7Tj' 6 events by maximizing the combined likelihood C. The statistical 
significance is estimated to be 1.7a by a comparison of the likelihood values for the best fit 
and a fit with zero signal yield. The E EC l distribution with the fit results superimposed 
is shown in Figure 03 The M(Ktt) distributions for the events in the signal region are 
also shown in Figure EJ while the momentum distributions of the reconstructed K*° — > Ktt 
candidates are given in Figure 03 




15 

miss 



5 -4 -3 -2 
Fisher Variable 



FIG. 2: The distributions for a) the momentum of K* candidate in the CM frame (P^-«), b) the 
missing momentum in the CM frame (-P^iss)' c ) ^ ne s< l uare d missing mass (M£ iss ), and d) the Fisher 
discriminant (arrow indicates the cut point). The dashed (solid) lines illustrate the distributions 
from signal (background) MC samples. 



TABLE I: The number of expected events in the signal and sideband regions. 



Signal region Sideband region 



Continuum e + e~ — > qq(q = u, d, c, s) 


0.5 ±0.5 


4.8 ±1.6 


Generic b — > c 


3.7 ± 1.4 


13 ±3 


Rare B decays 


0.5 ±0.2 


1.1 ±0.3 


Sum 


4.8 ± 1.5 


19 ±3 


Expected signal yield 


0.63 ± 0.03 


0.10 ±0.01 



Data 13 19 
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FIG. 3: The -Eecl distribution with fit results superimposed. The solid histogram illustrates the 
background distributions from MC simulations. The dashed, dotted and solid curves are the signal 
and background components and their sum, respectively. 

The systematic uncertainty assigned for MC statistics is 3.8%. By varying the PDF 
parameters and repeating the fit and observing the variation of the results, we infer an 
error associated with the fitting procedure of 5%. The charged kaon (pion) identification 
is checked by a high statistics D* tagged sample, and an error of 0.8% (1.2%) is assigned. 
Other possible systematic uncertainties such as tracking (2.1%) and number of BB (1%) are 
also included. The systematic uncertainty related to the K* mass selection is estimated to 
be 4%. The total systematic uncertainty is calculated to be 7.9%. Considering the effects of 
both statistical and systematic errors using an extension of the Feldman Cousins method , 
we obtain an upper limit of B(B° — > K*°uV) < 3.4 x 10~ 4 at the 90% confidence level. 

In conclusion, we have performed a search for the B° — > K*°vV decays with a fully recon- 
structed B tagging method on a data sample of 535 x 10 6 BB pairs collected at the T(AS) 
resonance with the Belle detector. We find 4.7lfg signal events with a statistical signifi- 
cance of 1.7 standard deviations. The obtained upper limit is three times more stringent 
than the only published constraint 0. Although the observed signal is much larger than 
the SM expectation, which is estimated to be 0.63 events in the signal region, the error is 
large, preventing us from drawing any conclusions at this point. We have examined several 
signal-like candidates and find that one of the events is consistent in missing mass with a 
B° — > K* 0/ y decay. The hard photon in that event hits the gap between the barrel and the 
forward endcap calorimeter. Further understanding of the events in the signal region will 
require much larger b — > c MC samples in addition to more data. The limit on B° — > K*°vV 
reported here is still one order of magnitude above the prediction of Buchalla et al. fH and 
hence still allows room for substantial non-SM contributions. 

We thank the KEKB group for the excellent operation of the accelerator, the KEK cryo- 
genics group for the efficient operation of the solenoid, and the KEK computer group and 
the National Institute of Informatics for valuable computing and Super-SINET network sup- 
port. We acknowledge support from the Ministry of Education, Culture, Sports, Science, 
and Technology of Japan and the Japan Society for the Promotion of Science; the Australian 
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FIG. 4: The M{Ktt) distributions from the events in the signal region, where a), b), c), d), and 
e) show the events from generic 6 — s- c MC, rare B MC, continuum e + e~ — > qq (q = u, d, c, s) 
MC, background MC combined, and data, respectively. The dashed histogram in d) illustrates the 
distribution for B° — > K*vV signal MC. 
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FIG. 5: The P*(Ktt) distributions (without Fisher requirement) from the events in the signal 
region, where a), b), c), d), and e) show the events from generic b — > c MC, rare B MC, continuum 
e + e _ — > qq (q = u,d,c,s) MC, background MC combined, and data, respectively. The dashed 
histogram in d) illustrates the distribution for B° — > K*vv signal MC. 
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